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Abstract 

We compute the dimension 6 effective Lagrangian arising from the tree level integration of an 

arbitrary number of bulk fermions in models with warped extra dimensions. The coefficients of the 

effective operators are written in terms of simple integrals of the metric and are valid for arbitrary 

warp factors, with or without an infrared brane, and for a general Higgs profile. All relevant tree 

level fermion effects in electroweak and flavor observables can be computed using this effective 

i , Lagrangian. 
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I. INTRODUCTION 

Models with warped extra dimensions [ll] offer a calculable path to the study of elec- 
troweak symmetry breaking (EWSB) induced by a strongly coupled sector j^j. Even if 
EWSB does not proceed precisely through the dual of the original Randall-Sundrum model, 
it is plausible that models with arbitrary warp factors provide us with a general enough 
sampling of the parameter space of strong EWSB models to give us the insight needed to 
decifer the LHC data. It is therefore crucial to develop techniques that allow us to study 
the phenomenological implications of models with warped extra dimensions for arbitrary 
backgrounds. 

Indirect constraints on models with warped extra dimensions have been extensively stud- 
ied in the past (see for instance [3|-l5|). Most of these studies are performed using a fixed 
background, and the few cases in which the low energy effective Lagrangian is computed 
for an arbitrary background, only the contribution of bulk gauge bosons has been taken 
into account. Fermion contributions on the other hand have been computed only for a 
fixed background, typically for a finite number of KK modes and always for specific models 
(using general results from intregration of four- dimensional vector-like fermions {(j, Q] , see 
also |8j). The purpose of this work is to fill this gap by computing the dimension 6 effec- 
tive Lagrangian generated after the integration at tree level of an arbitrary number of bulk 
fermions, including the effect of the full tower. The result is written in terms of integrals of 
the metric and the Higgs profile, thus being applicable to models with arbitrary warp factor 
and Higgs localization. EWSB effects can be included to all orders essentially only when 
they enter directly or effectively as boundary terms. In order to consider a more general 
case, we assume a light Standard Model (SM)-like Higgs with an arbitrary profile. This is a 
very good approximation in general if there is a light Higgs in the spectrum. For instance in 
composite Higgs models, bounds on the S parameter typically imply suppressed non-linear 
Higgs terms 4]. In order to obtain the background independent results, we have used 5D 
propagators to integrate out the full tower of bulk fermions (these techniques can be also 
useful for 5D loop calculations, see for instance Holographic methods would give the 



same results (see 



10[ for a detailed comparision of both methods). 



The rest of the article is organized as follows. We introduce our notation and discuss 
how to integrate out bulk fermions using 5D propagators in Section [Hj The different terms 
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entering the coefficients of the effective Lagrangian are computed in terms of integrals of the 
metric and Higgs profiles in Section III II An application of the formalism developed in these 
two sections is discussed in Section IIV} in which we compute the corrections to SM Yukawa 
couplings induced by bulk fermions. We conclude in Section |V] Technical details on the 
calculation and results of the 5D propagators and the general expression for the SM gauge 
and Yukawa couplings derived from the effective Lagrangian are given in two appendices. 



II. INTEGRATING OUT BULK FERMIONS 



In this section we will compute the effective Lagrangian resulting from the tree level 
integration of an arbitrary number of bulk fermions in models with warped extra dimensions. 
The only relevant assumption we make is to have a light Higgs in the spectrum, responsible 
for electroweak symmetry breaking, but otherwise remain completely general. Our methods 
can be generalized to include the effects of EWSB to all orders by including the Higgs 
vev as part of our general background or through boundary conditions (this is essential for 
instance in Higgsless models For clarity, we prefer to consider EWSB perturbatively 

as non-linear effects in the Higgs vev do not interfere with the effects of bulk fermions that 
we are interested in at the level of dimension 6 operators. Also, we will focus in this work 
on the quark sector but the results presented here are applicable, with minimal changes to 
the leptonic sector too (see [7|). With our assumption on EWSB, all fields can be classified 
according to their SM quantum numbers. In this case, the only new quarks that contribute 
to the effective Lagrangian at tree level and dimension 6 are the ones with the quantum 
numbers shown in table [H 



)(m) 
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TABLE I: New quark multiplets, Q {m \ that can mix with the SM fermions through Yukawa 
couplings. The index m labels the differen types of fermion multiplet additions in the given order. 
The electric charge is the sum of the third component of isospin T3 and the hypercharge Y . 
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We consider a 5D space with a general warped metric 

ds 2 = a 2 (z) [rj^dx 11 dx v — dz 2 ~\ , (1) 

where x M denote the standard four extended space-time dimensions and z parameterizes the 
extra dimension, which has two boundaries, L < z < L\, called the UV and IR branes, 
respectively. The warp factor is general but normalized to a(L Q ) = 1. Popular choices for 
this warp factor are a(z) = 1 for flat space and a(z) = Lq/z for AdSs. Also, the case of no 
IR brane can be simply obtained by taking L\ — > oo. 

Let us consider a number of bulk fermions, ^ b (x,z), classified according to their SM 
quantum numbers. Bulk fermions with quantum numbers other than the ones in table [I] 
do not contribute to the effective Lagrangian we are computing and are disregarded. We 
separate them into their zero mode (when it exists) and non-zero mode components 

*U*, *) = £(*) + *). ( 2 ) 

where, from now on, repeated indices imply summation unless otherwise stated, % — 1,2,3 
denotes the SM fermion generation index, and we have explicitly written the (4D) chiralities 

y L ,R = V L , R y, V l ,r=^^. (3) 

The action can then be written 

S 5 = S SM + AS, (4) 

where the part of the action involving only zero modes gives, after integration over the extra 
dimension, the SM action Ssm, for which we take for following convention, 

S SM = j d 4 x{qli0ql + u l R i0u l R + tf R i0d R - (V^\)<t L ju? R + \fql<Pd l R + h.c.) }, (5) 

whereas the part of the action involving heavy modes can be written 

AS = J d 4 x J 1 dz j (o b K 5 bc + 0^$ 6c ) ^ c + p J b + h.c] } , (6) 

where Ok represents the kinetic (including the SM gauge bosons through the covariant 
derivative) and mass contributions and O v the couplings to the SM Higgs. $ bc represents 
here the appropriate expression of the SM Higgs, according to the quantum numbers of the 
fermions involved. Explicitly we have 



[2y][l y+4 ] $ = 0, [3y][2 y+ x] -> $ = 0t_ 



[ 2y ][l y _x] -+ $ ^ 0, [3 Y ][2 Y _ k } $ = tf—, (7) 



where we have denoted in squared brackets the quantum numbers of the two fermionic fields 
involved in the coupling, cf) = i<j2<P* denotes the Y — — 1/2 Higgs doublet and the Pauli 
matrices are in the basis of T 3 L = ±1,0, the corresponding quark triplets. Finally, the 

currents J, which contain a SM fermion and the SM Higgs, have the generic form 

J b = ji + 4 = [ffM )s + f?M )s ] • (8) 

Due to chirality of the SM spectrum, only one of the two components is non-vanishing for 
each type of heavy fermion. In particular, using the notation in Table (TJ the non- vanishing 
currents can be written as 

Jl = -^(zWaftql JE = -X Dql (z)^ql (9) 

Jr = -\& r {z)4>Ur- (10) 

Jr = -^QXu' R (z)(pU R , jf = -X Q Y dR (z)4>d R , (11) 

JT = -K^W^^ Jt = -^(zWj—M, (12) 

where the X^^(z) encode the extra dimensional dependence of the effective Yukawa couplings. 
All these operators are functions of z that encode the dependence on the metric and the 
wave functions of the different fields. Their explicit expressions will be given in the next 
section. Their mass dimensions are [Ok] = 1, [O^] =0, [J] = 3 (also [^] = 2). 

We can integrate out the heavy fields at tree level by solving their classical equations of 
motion 

o b K s bc + e>J, c $H ^ c = - j b , (13) 

and inserting the solution back in the action 

AS = Jd 4 xdzJ b (x 1 z)^ b (x,z) = J -^- A dzJ a {p,z)¥{p )Z ), (14) 



where in the last equality we have switched to mixed position/momentum space [12| by 
Fourier transforming with respect to the four extended dimensions. We have implicitly 
denoted the functions and their Fourier transforms by their argument. \l/ b is supposed to be 
replaced with the solution of the equation of motion which, in mixed position/momentum 
space reads 

O b K (p, z)¥(p, z) = -J b (p, z)-j ^Of(z)<t> bc ( Pl )&(p - Pl , z). (15) 



Strictly speaking, one has to include the terms in Ok containing the SM gauge bosons as a 
convolution on the right hand side. Gauge invariance guarantees however that we can forget 
about those terms and recover them at the end by turning normal derivatives to covariant 
ones. We have explicitly checked that both methods give the same result. Let us define the 
Green's function of Ok with the zero mode subtracted 

O b K (p, z)P bc (z, z>) = 5 bc 5(z - z'), P bc = P bc - P b : TO mode , (16) 

which can be decomposed in their chiral (from the 4D point of view) components [3] 

Pll = V l PV r , P L r = V l PV l , P rr = V r PV l , P rl = V r PV r . (17) 

The different components can be written as follows 

P£(p; z, z') = (pL^PlKp; z, z'), P a R b R (p; z, z') = (fiLjP&fa z, z>), 

Pl b R { P - z, z') = P a L b R (p; z, z% P a R b L (p; z, z') = P ab L {p- z, z'), (18) 

with P(p; z,z') = (z, z') + {p 2 L\)P^ (z, z') + . . ., where we have introduced the appropri- 
ate powers of L\ to keep all components of the Green's functions with the same (vanishing) 
mass dimension. The solution of the equation of motion can be written in terms of these 
Green's functions by iteration 

**(p,z) = - J dz>P b %z,z'){j c (p,z>) + I ^0^(^)$ cd (p 1 )^(p-p 1 ,z , )} (19) 
= - J dz'P bc (z,z'){j c (p,z>) 

~ J dz" j ^<^(^(pi)P^ p1 (^^)^(p-Pi^' , )} + -- - • 
Inserting this solution back in the action and taking into account the following properties 

P hc * o =► j b r R = J b R Ji = o, (20) 
of ? o => J b L Ji = J b R j R = o, (21) 

we get the following effective Lagrangian 

£ 6 = -Li J dzdz' | J b L (x, z)P R °l bc (z, z')i0J c L (x, z') + (L i?)} 

+ I dzdz'dz"{j b L (x,z)P R f c (z,z')Of(z')^ (22) 



where we have switched back to position space (including the terms coming from the SM 
gauge fields), factored out the integration over 4D space to obtain the effective Lagrangian 
and the terms that we have not written are either absent due to the chirality of the SM 
spectrum or give rise to operators of dimension higher than 6. 

In order to write our effective Lagrangian in the standard basis of [14| we need to manip- 
ulate the operators in Eq. (|22|) . After applying the Leibniz rule for the covariant derivative, 
using the equations of motion from C,sm for the SM fermions and Fierz reordering, we end 
up with the following effective Lagrangian 6] 



4* = ^SM+^fE^Sj^^l^ + E^)^^ 



where we have explicitely factored out two powers of the cut-off A ~ L 1 1 (so that the 
coefficients of the different operators are dimensionless) , the sums run over all left-handed 
(LH) SM fields ipLi all SM fields ip and all SM right-handed (RH) fields ip R , respectively 
The different operators are defined as follows 

(Oflr = (0V^0)(^aVVi), (24) 

{O^f = (^ZVX^W), (25) 

{0 H f = (0 T e ^0)(^4), (26) 

(<D M ) ij = {tf4>){$L**Mi (27) 

where = <ft and &d R — 0- The first three kinds of operators come from the first line in 
Eq. (1221) and involve contributions from one kind of multiplet in table [I] at a time. The last 
operator, on the other hand, receives corrections from both lines in Eq. ( |22l) and therefore 
involves one or two different heavy multiplets at a time. The resulting coefficients of the 
effective operators are collected in tables HT1 and IHH written in terms of the following overlap 
integrals 

= -Li J dzd^X^^P^MX^^), (28) 
The only remaining pieces to compute the effective Lagrangian are the explicit expressions 



of the Otp(z) and X(z) functions and the zero momentum propagators. The former two are 
done in the next section whereas the latter are given in Appendix [A] 



TABLE II: Coefficients a™ resulting from the integration of an arbitrary number of each type of 
bulk quark. The coefficients are defined in Eq, ([28j) . 
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TABLE III: Coefficients a™ n resulting from mixing between bulk multiplets. The coefficients 

^ arc rlofintirl in TT.n 



7h*,j XX are defined in Eq. (|29j) . 
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III. GENERAL EXPRESSION FOR THE EFFECTIVE COEFFICIENTS 



Let us now give the explicit form of the different terms entering the coefficients of the 
effective Lagrangian. The relevant part of the action for the bulk fermions is given by 



S 



d A xdz a 4 ty b 



10+ (d z + 2-j 1 5 -aM^ b 



5 bc -af h \f^ hc \m c + S hoviIlA . 



(30) 



where Sbound. denotes any possible boundary term, b, c run over all bulk multiplets with 
quantum numbers appearing in table [J D denotes the SM covariant derivative, we have 
assumed a bulk Higgs with a zero mode profile given by fh( z ), normalized to 



3 f2 



dz a fi 



(31) 



Lo 



denote the 5D Yukawa couplings and $ bc represents the correct combination of the 
SM Higgs appropriate for the quantum numbers of ^b,c as described in Eq. (JTj). In the 
expressions above, we have assumed the Higgs to be a 5D scalar. In models of Gauge-Higgs 
unification, the Higgs comes from a gauge boson and the replacement fh — > or 1 fh should be 
made. Subtleties related to a boundary Higgs will be discussed in section IIV1 This action 
can be simplified with the following field redefinition 



where we have defined the effective (mass dependent) metric 

dM b (z) = exp — 2 / dz' a{z')M^ b {z' 
L Jl 

After this field redefinition, the action can be written 

S /(/'.a/:*" (i0alf b + 1 b d z y bc -aail 2 f h \f^ bc a 



(32) 



a -M b { Z >- 



(33) 



-7 5 /2 



+ Abound, (34) 



We separate now the bulk fermions into their zero mode and non-zero mode components 
as in Eq.© 

z) = fSSi^x) + *l R (x, z). (35) 

These zero modes realize the SM fermionic spectrum. We have denoted with i the flavor of 
the canonically normalized fermion zero modes. Note that we have allowed these fermion zero 
modes to be shared by several different bulk fields, as happens when boundary conditions 
mix different bulk fields. In the basis after the redefinition ( 132|) the zero mode profiles are 
just constants fixed by the boundary conditions and the orthonormality condition 



Qt 



dz a Mr 



n = f 



i 

R 



dz a M 



f u = f i 



dz a M r 



f D =5- 

a R 



(36) 



where as usual, a sum over all possible values of Q, U or D is understood. We can further 
require that the zero mode profiles diagonalize the zero mode Yukawa couplings for the down 



sector and the right part of the Yukawa couplings for the up sector. This amounts to 



Qt 



J 1 i L 



Qt 



2 2 



u 

J 

R 



f. 



D 

4 



(37) 



Note that for every zero mode profile that satisfy Eq. (I36I) we can define new ones 

s-Q . f Q T/t f [/ . fU xruR f D t-D irdR 



(38) 



with V L , V uR , V 3x3 unitary matrices so that the new profiles satisfy both (|36|) and (137} . 
This separation allows us to write the action in the form of Eq. fll]), where Eqs. ( 136H37]) 
guarantee that the SM action satisfies the convention in Eq. (j5J). The quadratic operator 
reads 



O 



'0ajf +j 5 d z 



+ bound, terms, 



(39) 



from which we can compute the relevant propagators as discussed in Appendix |A] The Higgs 
operator is given by 



O b ^ = -aa]^[ 2 f h \fJaJ 12 + bound, terms. 



Finally, the explicit expressions for the non-vanishing currents read 



(40) 





~ a f h(1 Mu X UQ a M Q f3 V jV 
Q L 


(41) 




= ^ afha M 2 n ^DQ a Mn f\ > 
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X Qd R { Z ) 
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X Q x u R { Z ) 


= Y^ a fh a M QX ^ < Qx u a M 2 u f^ 
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(45) 


^QYd R ( Z ) 
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(46) 




= a f h(1 M T x X T*Q a M Q f3 V jV 
Q 


(47) 


^TY q > L ( z ) 


= Y2 a fh a M T y X TYQ a M Q fj L V jv 


(48) 



Q 
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where possible contributions to boundary Yukawa couplings have not been explicitly written. 

We can now plug these explicit expressions, Eqs.f HOMHj) and (IA14HA29j) in Eqs. (I28|29p 
to compute the coefficients of the general effective Lagrangian for an arbitrary warp factor 
and Higgs profile. As a example of the use of these general equations we will compute in 
the next section the effect of bulk fermions in SM quark Yukawa couplings. 



IV. APPLICATION: HIGGS COUPLINGS 



As a straight-forward application of our formalism we compute flavour violating Yukawa 
couplings in the down sector of a simplified model with two bulk fermions with the quantum 
numbers of Q and D (following the notation in Table H]) and boundary conditions Q ~ 

[++], D ~ [ ]. The two signs denote the boundary condition at the UV and IR brane, 

respectively, and a + (— ) denotes Dirichlet boundary conditions for the RH (LH) component 
of the bulk field at the corresponding brane. Their zero modes give rise to the SM 
and da quarks. We want to compute the modified SM Yukawa couplings that result after 



the integration of the heavy fields. This has been studied recently in Refs. [15|, [l6j. In 
particular, Ref. clarifies the implications of brane-localized Yukawa couplings for fields 
with Dirichlet boundary conditions. As we will see, the use of 5D propagators makes it 
apparent the ambiguity of these couplings and a very simple regularization of the brane 
terms gives directly the same result as previously obtained with other methods. 1 The SM 
fermion gauge and Yukawa couplings derived from the effective Lagrangian in Eq. (123]) were 



computed in [6J. We collect the main results in Appendix |BJ The modified Yukawa couplings 
can be written, in the physical basis {i.e. with diagonal fermion masses including effects of 
order v 2 /A 2 ) 

C H = --^(C^X + d\Y^ R )E + h.c, (49) 



1 Similar ambiguities occur in brane localized kinetic terms for fields with Dirichlet boundary conditions 17 1. 
In that case, there is a well-defined prescription to deal with such ambiguities in terms of field redefinitions 
and classical renormalization [18fl . 
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where the coefficients read 

,2 



VV; = <),,X", - — - I Vik(atu<j>)kj + -78ij[Vik(ot u <i>)kj + { a u<t>)\ k VL 



v 

V 2 ( 1 



(50) 



with v « 246 GeV the Higgs vacuum expectation value. Note that we have not written 
couplings proportional to d^H because they vanish in our case due to the hermiticity of 
the corresponding operator coefficients (see Appendix [Bj. These modified couplings are thus 
fully determined by two coefficients that read, in our example, 



A 2 2 



1 R Ul R 



1 nD;R i 1 \dnQ;L , DQ;RL 

A 2 = 2^ + 2 K ^ R < + \< ■ (52) 

Using the explicit expressions for the propagators in Appendix we can compute the 
corrections to the Yukawa couplings for arbitrary backgrounds and Higgs profiles. In the 
case of an IR boundary Higgs, the terms proportional to /3 contain only Yukawa couplings 
of fields with Neumann boundary conditions. They are well defined in the thin brane limit 
and a direct application of our formulae reproduces the results in the literature. The term 
proportional to 7 on the other hand, involves a Yukawa coupling of fields with Dirichlet 
boundary conditions and therefore should a priori vanish. We will show that the result is 
ambiguous in the thin brane limit but a simple regularization of the brane gives a result 
that agrees with the ones presented in the literature. 

In order to directly compare with the results in the literature we consider a pure AdSs 
background, a boundary localized Higgs and a constant bulk mass for the fermions M = 



c/Lq. Following 15( we define the Yukawa Lagrangian for a boundary Higgs (for simplicity 



we focus on the down sector) 

S Yuk = - / d 4 xdzS(z - Li)a 3 \Q L Y* D L Q (j>D R + Q R Y^ D L (j)D L + h.c. + . . . , (53) 



where we have already canonically normalized the Higgs boson (so that its vev is v ~ L^ 1 ) 
and we have included a term involving boundary values of fields with Dirichlet boundary 
conditions. The latter terms do not involve the zero modes and therefore do not affect the 
values of the effective z-dependent Yukawa couplings appearing in the currents, that read 
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m our case 



\ Dq {z) = 6(z - L^KY^y XioL/, 



A Qd (z) = Siz-LjalY^L^J^ 



cD 



However, both terms contribute in principle to the O b £ term, which reads in our case 



8(z - L X )L X a M 2 Y* D] a 2 M V L + a 2 Mo Y 2 01J1 a M 2 V 



5Df 



R 



(54) 
(55) 

(56) 



In fact, due to the chiral structure of the coefficient 7 I ^j' RL only the term proportional to 
Y 2 can give a non-vanishing contribution. Now that we have the explicit expression for the 
relevant overlaps with the Higgs, we only need the corresponding propagators to compute 
the modified Yukawa couplings. The /3 coefficients are given by 



IlIl 
ft* 



-L\f? L ] aiYl> D aJ D P^l'\c Dj L u ^a^Y^aljg, 

i 



(57) 

L lfdR a M D Y l D]a ll Q PLL ++ \ C Q^ L h L l) a M Q Y l Da M D fd R i ( 58 ) 



where aj/(ii) = (L / Li) 2ML ° , cq = MqL , cd = Me>L and the explicit expressions of the 
boundary propagators for the AdS case is 



P, 
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+ (1 



X 



(59) 



(3-2c)(| 



2c-3 



See Eqs. (IA14|) and (lA23j) . The 7 coefficient on the other hand, involves propagators that 
present discontinuities at the branes 



P, 



(0)[H 



RL 



\z,z') = 6{z-z')- 



dzidMizi, 



and 



1 



z -z 



dz±a-M(zi). 



(60) 



(61) 



This discontinuities indicate that the coefficients 7 are ambiguous in the thin brane limit 
and therefore the Dirichlet brane Yukawa term, proportional to Y 2 , can potentially give a 
non- vanishing contribution. In order to obtain a finite result we will replace the brane terms 
in Eqs. ( |54ll56l) with regulated delta functions 



S(z - L\) -> 5 e (z - Li) 



z < Li — e, 

Li — e < z < L\. 



(62) 



13 



Once the branes are regulated, we can perform all the relevant integrals and then take the 
limit e — > 0. The (non- vanishing) result we obtain is 

7 Sf = ^f?>k Y " DY ^ Y " Da tfl- (63) 
Using the fact that (see Appendix 



- j0 = ^<^, (64) 

the result for (at^)y deduced from Eqs. fl52|63p exactly reproduces the results of Eq. (73) in 
Ref. 15], taking into account the following dictionary between our notation and theirs: 

v h - = V2v^, (ai Q fgr™ = (F q )W, (a'Jj d D R ) herc =(F d f 5] . (65) 

In Ref. [l5| the correct result was obtained, either by taking the boundary limit of a bulk 
Higgs or by summing the effect of a number of KK modes of the order of the inverse width 
of the regulated brane in the case of a brane Higgs. In our case, this resummation is done 
automatically by means of the 5D propagators which immediately shows the ambiguity 
present in terms of Dirichlet Yukawa coupligns and the need for a regularization of such 
couplings. 



V. CONCLUSIONS 



Models with warped extra dimensions represent calculable examples of models in which 
EWSB proceeds through a strongly coupled conformal sector. With the LHC currently 
probing the TeV scale, it is important to have a handy way of computing the low energy 
effects of the largest possible number of different models of strong EWSB. The effective 
Lagrangian of models with warped extra dimensions arising from the tree level integration 
of bulk gauge bosons has been know for general background metrics and Higgs profiles for 
some time now. In this work we have completed the most relevant part of the tree level 
(dimension 6) effective Lagrangian by including the effect of an arbitrary number of bulk 
fermions, again with general warp factor and Higgs profile. Our general calculation involves 
the integration of the full 5D fields (with the zero modes subtracted when present) by means 
of the 5D propagators. This has the advantage of trading eigenvalue searches and sums with 
integrals that can be easily done numerically. Also, it makes apparent the subtleties related 

14 



to brane localized terms involving fields with Dirichlet boundary conditions as the automatic 
resummation of all the modes gives a non-vanishing propagator near the brane for certain 
components of fields with Dirichlet boundary conditions. 

Our results are completely general and can be applied to any model with warped (or 
flat) extra dimensions, with or without an IR brane 19J. This finally permits the complete 
calculation of the low energy effects of a large variety of holographic models of strong EWSB, 
including all the relevant electroweak and flavor effects. 



Note Added: 



Upon completion of this work, Ref. [20| appeared in the arXives. In that reference, the 
contribution of bulk fermions to the coupling of the down sector to the Z boson is computed 
using similar techniques to the ones presented here. 



Acknowledgments 

We would like to thank Manuel Perez- Victoria for useful discussions. This work has been 
funded by MICINN (FPA2006-05294, FPA2010-17915 and FPU programs) and Junta de 
Andalucia (FQM 101, FQM 03048 and FQM-6552). 

Appendix A: Derivation of the fermionic Green's functions 

Assuming no boundary terms, the quadratic operator in mixed momentum/position space 
reads 

(Al) 

The different components of the fermionic Green's functions satisfy the following coupled 
equations 

al}P RL - L x d z P LL = 0, a M L lP 2 P LL + d z P RL = 5{z - z'), (A2) 
a M P M + ii^ = 0, aj}L lP 2 P RR -d z P LR = 5(z-z'), (A3) 



O 



K 



'■M 



+ l 5 d z 
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and can be written in terms of the KK mode profiles as follows 

/^)/ n V) 



9 9 ? 

p z — m z n 

9 9 ' 

p z — m~ 



p 2 — 



(A4) 
(A5) 
(A6) 



We can obtain equations for the LL and RR components by iteration 



Li 



d z a M d z + a M p 2 Pp L (z } z') = 5(z - z'), 



Ll 



d z a. M d z + a. M p 2 P« R (z,z') = S(z-z'). 



(A7) 
(A8) 



These equations are identical to the ones of a gauge boson with a generalized metric cl±m(z) 
for LH and RH components, respectively. The zero momentum propagators for gauge bosons 
with an arbitrary warp factor can be computed in a number of ways. The most direct method 
is probably the one suggested recently in which amounts to solving directly the equations 
that the Green's function at zero momentum satisfy. The resulting Green's functions at zero 
momentum, with zero modes subtracted and for an arbitrary warp factor a, read (see 



1 i 

L 



1 



Li 



+ 



L 2 , L 



dz 2 a~ l {z 2 ) / dzxa(zx) 
l Jl L 

Li pzi pLi 

dz\a~ {zi) / dz 2 a(z 2 ) / dz 3 a{z 3 

o J Lq J Z\ 



dz 2 a 1 (z 2 ) 



LiPWt— ](a,z, z') 

L 1 P<M- + \a,z,z') 
IaPW+^faz, z') 



f^dzta l {z x ) f^dz 2 a 1 (z 2 ) 



z< 

Lo 
Li 

z> 



f L *dz 3 a \z 3 ) 
dz\a~ (zi), 

dz\a~ (zi) , 



dz\a\Z\) 
(A9) 

(A10) 
(All) 
(A12) 



where the superscript denotes the boundary conditions with a + (— ) denoting Neumann 
(Dirichlet) boundary conditions and the first (second) sign refering to the UV (IR) boundary 
and we have defined 

/•Li 

L= dza(z). (A13) 

Jl 

The fermionic Green's functions P^l and P^ can be expressed in terms of these with the 
replacement a — > a±M (and the understanding that the boundary conditions correspond to 
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the LH chirality of the 5D field). The LR and RL components can then be trivially obtained 
from these by using the first equations in (IA2I) or (IA3I) . We can also use Plr{z,z') = 
Prl(z',z). 

Let us now give the result for the different possibilities of boundary conditions. 

- [++] case 

The LL and RR components are given by 

Pi%,z') = P^(a M ,z,z'), (A14) 
P { R %(z,z') = P^-\a„ M ,z,z>), (A15) 

whereas the mixed ones read 

Pj®(z, z') = 6(z' - Z )--L [ dzxauizx), (A16) 

L M JLo 
1 f* 

z ) - — / dz x a M {z x ), (A17) 



P£l(z,z') = 9(z 
[ — h] case 



P?l(z,z>) = P^-+\a M ,z,z'), (A18) 

P R %(z,z') = P^ + -\a_ M ,z,z'), (A19) 

P^(z,z') = -e(z-z'), (A20) 

P R °l(z,z') = -e(z'-z), (A21) 



case 



Pl°l(z,z') = PW-\a M ,z,z'), (A22) 

P R %,z') = P^ + +\a„ M ,z,z>), (A23) 

P®(z, z') = -9(z - z') + J- f dz 1 a~ M 1 (z 1 ), (A24) 

P R % } z') = -0(zf - Z ) + J— [' dz 1 a^(z 1 ), (A25) 

L -M Jl 
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- [+— ] case 



P^(z,z') = P^ + -\a M ,z,z'), (A26) 

P R %(z,z') = P<M-*(a-M,z } zf), (A27) 

PjgL(z,z>) = e(z>-z) } (A28) 

P^(z,z') = e(z-z'). (A29) 

Appendix B: Trilinear couplings in the physical basis 

In this appendix we reproduce the couplings of the SM quarks to the SM gauge bosons 

and Higgs with the Lagrangian in Eq. (|23|) that were computed in [6]. The masses of the 
SM quarks in the physical basis read 



m "(phys.) 



-j= f A? - -[(Va^u + (Vo^)],]^ , (Bl) 



? phys ° = 7f ( A " - + MU^) ■ ( B2 ) 



In this basis, the SM quark couplings to the SM gauge and Higgs bosons can be written as 



-2 L X*"fj L - J^f *fd? R - 2 sin 2 6 W J£ U ) Z m 
C w = -A={uiWt^d j L + u R W^d j R )W+ + h.c, (B3) 

+ 7X + u R Zf^u R - dlZf^di - d R Zf^d R ) id,H. 
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The unbroken U(1)q protects the terms proportional to J£ M . The expressions to order 1/A 2 
of the coupling matrices X, W, Y and Z in terms of the coefficients a x are: 

A ij ~ d H~ Vik ^ a <Pi +a H ~ a <t>q ~ a H > klV ij> 
1 v 2 

^ = ^ + ^(«2+«2 t +«2+«S? t )«. 

x f = ( B4 ) 

= Vik (Ski + ^(oQ)ki 
1 v 2 

n 



ij ~~ _ 2A2 ^ ^ ~ a H ~ a H +a H )kl V lji 
Z dL - -—(a {1) - a (1)t + a (3) - a (3)t V- 

We have introduced the unitary matrix 



v 2 



with 



V = V + -(VA d L -AlV), (B5) 



1 LjiJ 2 l 2^' J (Af)2_(_i)^(A^) 2 ' 

(B6) 
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Note that, to order 1/A 2 , we can substitute V by V everywhere in Eq. (1B4I) . so that the 
different couplings depend on only one unitary matrix. 
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